Coherently induced grating in refractive index enhanced medium
Liu Zhuan-Zhuan, Chen Yu-Yuan, Yuan Jia-Yu, Wan Ren-Gang
School of Physics and Information Technology, Shaanxi Normal University, Xi’an 710062, China

 

† Corresponding author. E-mail: wrg@snnu.edu.cn

Project supported by the National Natural Science Foundation of China (Grant Nos. 11204367 and 61475191) and the Fundamental Research Funds for the Central Universities, China (Grant Nos. GK201503022 and GK201601008).

Abstract

We demonstrate a scheme for coherently induced grating based on a mixture of two three-level atomic species interacting with two standing-wave fields. As a result of interaction between the absorptive and amplified Raman resonances, the refractive index of the medium can be enhanced and modulated periodically. Then a sinusoidal grating, which can diffract the probe field into high-order directions, is coherently formed in the medium. The proposed scheme is theoretically investigated in a mixture of atomic isotopes of rubidium. The results show that the diffraction efficiency depends strongly on the two two-photon detunings of the two Raman transitions and the intensities of the two driving standing-wave fields. The proposed electromagnetically induced grating scheme may be applied to the all-optical switching and beam splitting in optical networking and communication.

1. Introduction

Recently, the research of light-media interaction has played a significant role in the area of quantum optics, where many interesting optical phenomena can be caused by atomic coherence and quantum interference. Of them, the electromagnetically induced transparency (EIT) is the most typical.[1,2] The EIT is mainly used to manipulate the refractive and dispersive properties by suppressing or eliminating the absorption for a weak probe beam in atomic medium driven by an additional strong coherent field. Meanwhile, coherent population trapping (CPT) state,[3,4] or the so-called dark state, can be formed in this process. The principal mechanisms for EIT and CPT have been explored in many recent applications, such as lasing without inversion,[5,6] refractive index enhancement,[7,8] slow light,[9,10] light storage,[11,12] enhanced Kerr nonlinearity,[13,14] optical solitons,[15,16] etc.

As is well known, an electromagnetically induced grating (EIG) can be created in a medium by replacing the travelling-wave driving field in EIT with a standing-wave. The probe field will obtain a periodic variation of absorptive coefficient due to the spatial modulation induced by the standing-wave field. At the same time, the refractive index also changes periodically.[17] Then the medium is equivalent to a grating on which the resonant probe beam can be diffracted into high-order directions. This phenomenon has been observed experimentally in cold and hot atomic samples[18,19] and universally applied to optical bistability,[20] all optical switching and routing,[21] beam splitting and fanning,[22] quantum Talbot effect,[23,24] surface solitons,[25] photonic Floquet topological insulators,[26] etc. However, the diffraction efficiency of EIG is not enough due to the inevitable atomic absorption. Therefore, an electromagnetically induced phase grating (EIPG) was proposed for high diffraction efficiency by using enhanced Kerr nonlinear refractivity at low light level.[27,28] Likewise, one can build gratings in other media or configurations such as quantum well,[29] quantum dot,[30] molecular magnets,[31] double-dark resonances system,[32,33] phase-dependent system[34] and CPT system.[35] Meanwhile, a kind of gain-phase or gain grating based on the effect of active Raman gain (ARG) has been put forward to obtain much higher diffraction efficiency in high-order directions.[36] Recently, Wang et al. extended EIG to the two-dimensional (2D) case and achieved an electromagnetically induced cross-grating (EICG) in a four-level tripod-type atomic system.[37] Then, schemes for EICG based on EIT or ARG in multi-level atomic systems were proposed, and the diffraction efficiencies were improved.[3841]

In this paper, we theoretically investigate a one-dimensional (1D) phase grating coherently induced in a mixture of two species of Λ-type atomic system with EIT and ARG, respectively. This three-level atomic system has been proposed for the enhancement of refractivity without absorption.[42,43] Recently, the creation of a spatially distributed PT-symmetric refractive index was explored.[44] Here, by using two far-detuned standing-wave fields to drive the corresponding atomic transitions in the vicinity of two-photon resonance, we obtain spatial absorption-phase, pure-phase, and gain-phase gratings through modulating the interaction of the two Λ schemes. Thus, the probe beam can be diffracted to high-order directions with different intensities. It is found that the diffraction efficiency of the grating depends on the detunings of the two two-photon Raman transitions and the amplitudes of the standing-wave fields. We also analyze the physical mechanism for such a mixture of EIT and ARG in a real atomic system with 85Rb and 87Rb. The scheme can be used as all-optical switching and beam splitter which have potential applications in optical networking and communication.

2. Model and equations

The considered atomic system is a mixture of two three-level Λ-type atomic species (EIT and ARG) interacting with three coherent laser fields as illustrated in Fig. 1. The atomic densities of the species are N1 and N2, respectively. Then the total atomic density is N = N1 + N2. Each of Λ systems is marked by si (i = 1,2) with one excited state |3,si〉 and two ground levels |1,si〉 and |2,si〉. A weak probe field with frequency ωp and a pair of coupling fields with frequencies ωs1 and ωs2 drive the corresponding transitions. The Rabi frequencies of the driving fields are denoted as Ωsi = μlm,siEsi/2ħ and Ωpi = μlm,siEp/2ħ (i = 1,2), where Esi and Ep represent the amplitudes of electric fields and μlm,si (l, m = 1, 2, 3) are the dipole moments of the corresponding transitions. We define that Δs1 = ω31,s1ωs1 and Δs2 = ω32,s2ωs2 are the one-photon detunings for the two coupling fields, and Δp1 = ω32,s1ωp and Δp2 = ω31,s2ωp are the detunings for the probe field with ωlm,si(l,m = 1,2,3; i = 1,2) being the resonant frequencies of corresponding transition in si-system. Here, all the involved fields are detuned far from the excited state |3,si〉. For the EIT system, the atom is transparent to the probe field under the condition of two-photon resonance, i.e., δ1 = 0, and its refractive index is unity. However, one can obtain a narrow absorption peak near the two-photon resonance.[45] At the same time, the refractivity is enhanced due to the Raman process. For the ARG system, the probe field can be amplified by the atoms and a nonzero susceptibility is then achieved when the two-photon detuning δ2 ≠ 0. The combination of Raman absorption and Raman gain can result in enhanced refractivity with vanishing absorption. The laser driven 85Rb and 87Rb atoms can be utilized as the two three-level Λ systems, of which one provides an absorption resonance and the other induces a gain resonance for the probe beam.[42,43]

Fig. 1. (color online) Energy level scheme based on a mixture of two species of Λ atomic system driven by three light fields. The initially populated levels are indicated by the dots.

In the framework of the semiclassical theory, with the rotating-wave approximation and the electric-dipole approximation, the Hamiltonian of si-system in the interaction picture can be written as In Eq. (1), the quantities δ1 and δ2 are the two-photon detunings of the laser beams from each Raman transition, and they are defined as δ1 = ωs1ω21,s1ωp and δ2 = ωs2 + ω21,s2ωp, respectively.

By using the Liouville equation, we obtain the relevant density-matrix equations as follows: with and . In the above equations, we consider the relaxation process, and γlm,si represents the dephasing rates of |l,si〉 ↔ |m,si〉 transitions. We assume that the atoms are initially populated in one of the two ground-state sublevels, i.e., ρ22,si = 1(i = 1,2).

By solving Eqs. (2a)–(2f) in the steady state, the coherent terms corresponding to the probe transitions can be derived as follows:

In this system, the refractive index enhancement of the medium relies on the interaction strength between two Raman transitions, the first of which results in two-photon absorption and the second one contributes two-photon gain to the probe field. Therefore, the total complex susceptibility is equal to the sum of individual contributions from each of the Λ-type schemes, and the absorption and refractivity can be controlled by changing the two-photon detunings δ1 and δ2. According to the polarization of the medium, Pp = ε0χEp = N1μ23,s1ρ32,s1 + N2μ13,s2ρ31,s2, the susceptibility of the weak probe field can be written as where ε0 is dielectric constant in a vacuum.

Therefore, the susceptibility of the whole system is given by

In order to investigate the coherently induced grating in refractive index enhanced medium, we use two standing-wave fields Ωs1(x) = Ωs2(x) = Ωssin(πx/Λ) to drive the atoms, with Λ being the spatial period, and assume that the probe field propagates along the z direction through an atomic medium with length L. The propagation of probe field determined by the atomic polarization, which, in the slowly varying envelope approximation and the steady-state regime, can be described by the Maxwell’s equation as where λp is the wavelength of the probe field. By solving Eq. (6), the transmission function of the probe field at z = L can be derived as The terms of α(x) = πIm[χ]L/λp and ϕ(x) = πRe[χ]L/λp represent the absorption coefficient and phase shift of the probe field, respectively. According to the Fraunhofer diffraction theory, we obtain the Fraunhofer-diffraction intensity as where corresponds to the Fraunhofer diffraction of single spatial period, θ is the diffraction angle with respect to the z direction, and M is the number of spatial periods of the grating. Through Eq. (8) with grating equation sin θ = p/Λ, which determines the n-th order diffraction angle, we can calculate the zeroth-order, first-order and second-order diffraction intensities as follows:

In the following, we focus on a real system and consider two Λ schemes in a dual magneto-optical trap of rubidium atomic isotopes with 85Rb (for EIT) and 87Rb (for ARG) using the D1 line (λp = 795 nm), and assign |1,s1〉 = |52S1/2, F = 2〉, |2,s1〉 = |52S1/2, F = 3〉, and |3,s1〉 = |52P1/2, F = 3〉 for s1-system, and |1,s2〉 = |52S1/2, F = 1〉, |2,s2〉 = |52S1/2, F = 2〉, and |3,s2〉 = |52P1/2, F = 1〉 for s2-system. In the following calculations, the parameters of this atomic system are N1 = N2 = 5.0 × 1012 cm−3, μ23,s1 = μ13,s2 = 2.5377 × 10−29 C · m, γ = γ32,s1 = γ31,s2 = γ32,s2 = 5.75 MHz/2, and γ = γ12,s1 = γ12,s2 = 0.05γ. We can manipulate the susceptibility by changing the two-photon detunings δ1 and δ2 to obtain diverse phase gratings, and the variations of intensities of driving fields can be used for achieving different phase shifts to enhance the diffraction efficiency.

3. Results and discussion

In order to acquire a high diffraction efficiency for probe beam with the coherently induced grating, we must pay attention to two factors: a large phase shift should be firstly possessed across the probe beam, and the transmission of the probe must be large enough. By using the equations as outlined above, we investigate the formation of phase grating and its diffracting power in this section, and discuss the Fraunhofer diffraction patterns of the probe field in three types of gratings with different phase shifts.

The susceptibility of the medium plays a significant role in the formation of EIG, and we firstly focus on the probe absorption Im(χ) and dispersion Re(χ). As discussed above, we can efficiently modulate the strength of the two-photon Raman absorption and Raman gain via manipulating the two-photon detunings δ1 and δ2 of the laser beams, so that the absorption can be reduced and the probe can even be amplified. At the same time, the real part of the susceptibility also obtains a large value arising from the constructive contributions of two Raman transitions. In addition, the separation value δ1δ2 of two Raman resonances determines the interaction between EIT and ARG systems, and then the refractive index is dependent on the separation of Raman resonances.[42,43] For δ1δ2γ, two Raman transitions are independent of each other and Re(χ) drops. But for δ1δ2γ, the effect of two Raman transitions can mutually offset. In such a case, both the absorption and the refraction are cancelled. Therefore, we choose that the separation value of two Raman resonance satisfies δ1δ2 = 3γ here, which is of the order of γ.[42,43] As shown in Fig. 2, the dispersive part Re(χ) and absorptive part Im(χ) of the probe susceptibility are plotted separately as a function of two-photon detuning δ1 for δ1δ2 = 3γ, marked with the blue-solid line and red-dashed line, respectively. We take the Rabi frequencies of two driven fields as Ωs1 = Ωs2 = 4γ. We consider three different conditions for which δ1 = 0.5γ (point A), δ1 = 1.5γ (point B), and δ1 = 2.5γ (point C), respectively. At point B, the absorptive part of susceptibility [Im(χ)] vanishes due to the balance between two-photon absorption in EIT system and two-photon gain in ARG system while the dispersive part [Re(χ)] is enhanced due to their constructive contributions. For δ1 < 1.5γ, the atomic medium is absorptive to probe beam. In this case, the slight deviation from the EIT resonance results in an absorption which is stronger than the gain provided by ARG. However, under the condition of δ1 > 1.5γ, the Raman gain is dominant and the probe light can be amplified by the atoms. Nevertheless, Re(χ) is always enhanced, that is, the atomic system is a refractive index enhanced medium (see points A and C).

Fig. 2. (color online) Real (solid curve) and imaginary (dash curve) parts of the probe susceptibility versus the two-photon detuning δ1 under the condition of δ1δ2 = 3γ. Parameters are N1 = N2 = 5.0 × 1012 cm−3,γ = γ32,s1 = γ31,s2 = γ32,s2 = 5.75/2 MHz, γ = γ12,s1 = γ12,s2 = 0.05γ, Δs1 = −862γ, Δs2 = 862γ, and Ωs1 = Ωs2 = 4γ.

The strength of the two-photon absorption in EIT system and the two-photon gain in ARG system are dependent on the amplitudes of the coupling fields [see Eq. (5)]. Then coupling fields with standing-wave pattern can result in a spatially modulated grating for the probe fields. As shown in Fig. 3, the properties of grating are different for the probe field when two-photon detunings of the two Raman transitions δ1 and δ2 take different values. With the parameters fixed, a large phase modulation is observed, reaching a peak phase-shift value of π at x = ±0.5Λ for Ωs = 4.5γ [see the solid lines in Figs. 3(a1)3(a3)]. However, the amplitude modulations vary with different two-photon detunings [see the dashed lines in Figs. 3(a1)3(a3)]. Figure 3(a1) shows the transmission function for the two-photon detunings of δ1 = 0.5γ and δ2 = −2.5γ. There is a certain amplitude modulation of transmission function, which oscillates around an average transmissivity of 46.8% along the x direction in a period of Λ. Therefore, the probe field has a significant energy loss inside the atomic medium due to the absorption effect in the EIT system with non-zero two-photon detuning. This is a limit to the diffracting power of the atomic grating. The resulting first-order diffraction efficiency of grating reaches only to 8.9% as shown in Fig. 3(a2). Thus, this grating is a mixture of an absorption grating and a phase grating. Figure 3(b1) shows the transmission function for δ1 = 1.5γ and δ2 = −1.5γ. We can obtain a phase modulation of the order of π for the probe beam with extremely low energy loss in the atomic medium, which is due to a good balance between the two-photon absorption in EIT system and the two-photon Raman gain in ARG system. As a result, a pure phase grating with nearly 100% transmissivity is developed in the medium and the diffraction efficiency of probe beam is improved. The energy diffracted into the first-order direction can reach to 31.5%, which means that this EIG is very close to an ideal sinusoidal phase grating (about 34%) as shown in Fig. 3(b2). In the same way, we plot the transmission function with the detunings being δ1 = 2.5γ and δ2 = −0.5γ in Fig. 3(c1). It shows that the amplitude modulation of transmitted probe light is periodically modulated in space with an average transmissivity large than 100% which corresponds to an amplification of the probe beam. At the same time, the resulting phase modulation within a period is still π. Here, the probe gain is attributed to the effect of ARG. The far-field diffraction pattern of this grating is depicted in Fig. 3(c2). It is clear that the first-order diffraction efficiency can exceed 100%. Therefore, a kind of gain-phase grating is obtained. From Figs. 3(a2), 3(b2), and 3(c2), we can observe that a small amount of energy is also shifted into the second diffraction orders which are located at sin θ = ±0.5, and the diffraction efficiency is enhanced with the increase of total transmittance.

Fig. 3. (color online) Amplitude |T(x)| (red dashed line) and phase ϕ(x) (blue solid line) of the transmission function T(x) plotted over two space periods for different two-photon detunings: (a1) δ1 = 0.5γ, δ2 = −2.5γ; (b1) δ1 = 1.5γ, δ2 = −1.5γ; (c1) δ1 = 2.5γ, δ2 = −0.5γ. ((a2), (b2), and (c2)): Plots of diffraction pattern I(θ) versus sin θ corresponding to the (a1), (b1), and (c1) cases. Parameters are Ωs = 4.5γ, L = 3 mm, M = 6, and Λ = 4λp. Other parameters are the same as those in Fig. 2.

The diffraction efficiency depends on the two-photon detunings δ1 and δ2 and the intensity of driving standing-wave field Ωs. Figure 4 displays the changes of transmission function and diffraction pattern for Ωs = 5.8γ. Other parameters are the same as those in Fig. 3. It is obvious that we can obtain a phase modulation of the order of 1.5π at x = ±0.5Λ [see the solid lines of Figs. 4(a1), 4(b1), and 4(c1)]. However, the amplitude modulation |T(x)| varies with two-photon detuning. For the cases of δ1 = 0.5γ and δ2 = −2.5γ, the amplitude modulation is plotted in Fig. 4(a1) with the dashed line. Its average transmissivity is about 38.2%. Hence, the medium absorbs more energy from the probe beam than that of the case in Fig. 3(a1). It is due to the increased two-photon absorption in EIT system. The diffraction intensity corresponding to this absorption-phase grating is plotted in Fig. 4(a2). Comparing Fig. 4(a2) with Fig. 3(a2), we can see that more probe energy is deflected into the second order diffraction peaks. This is due to the improvement of phase modulation. When the detunings are δ1 = 1.5γ and δ2 = −1.5γ, we can see that there is no amplitude modulation as shown in Fig. 4(b1) with the dashed line. Consequently, a pure-phase grating is built up in the medium. Almost all the energy in the probe beam is diffracted into the first-order and second-order directions with nearly no light left in zeroth-order [see Fig. 4(b2)]. The diffraction efficiency is enhanced by comparing with Fig. 3(a2). In Figs. 4(c1) and 4(c2), we also plot, respectively, the transmission function |T(x)| and diffraction pattern I(θ) for detunings δ1 = 2.5γ and δ2 = −0.5γ. There is a periodical gain modulation for probe field, which leads to a structure of gain-phase grating, and therefore the diffraction efficiencies of high-order directions are significantly enhanced. As shown in Fig. 4(c2), the corresponding values of each high-order diffraction intensity increase to I(θ1) = 450% at sin θ1 = ±0.25 and I(θ2) = 306% at sin θ2 = ±0.5.

Fig. 4. (color online) Amplitude |T(x)| (red dashed line) and phase ϕ(x) (blue solid line) of the transmission function T(x) plotted over two space periods for different two-photon detunings: (a1) δ1 = 0.5γ, δ2 = −1.5γ; (b1) δ1 = 1.5γ, δ2 = −1.5γ; (c1) δ1 = 2.5γ, δ2 = −0.5γ when Ωs = 5.8γ. [(a2), (b2), and (c2)] Diffraction patterns I(θ) as a function of sin θ for the (a1), (b1), and (c1) cases. Other parameters are the same as those in Fig. 3.

Additionally, comparing the corresponding graphs in Figs. 3 and 4, it is found that the depth of phase modulation and the diffraction efficiency are dependent on the Rabi frequency of coupling field. The phase modulation makes the probe light mainly disperse into the high-order directions. Meanwhile, the amplification of probe beam from the Raman gain can enhance the diffraction efficiencies of all orders.

In order to further clearly clarify the effect of coupling field intensity on diffraction efficiency, we display the variations of the different order diffraction intensities I(θ0), I(θ1), and I(θ2) with coupling field intensity Ωs for different values of δ1 and δ2 in Fig. 5. It is found that the I(θ0) decreases as Ωs increases as shown with the solid lines in Figs. 5(a) and 5(b) for the absorption-phase gratings and phase gratings. But for the case of gain-phase grating (δ1 = 2.5γ and δ2 = −0.5γ, the zeroth-order diffraction intensity I(θ0) firstly presents two gain peak at a certain Ωs, beyond which it declines for large Ωs as shown with solid line in Fig. 5(c). Then the high-order diffraction intensities I(θ1) and I(θ2) can successively increase as Ωs further increases, and they can take their maximum values under an appropriate Ωs with the dashed and dashed-dotted lines shown in Figs. 5(a)5(c), respectively. Moreover, the maximum values of I(θ1) and I(θ2) are not equivalent to each other for different conditions of δ1 and δ2. As a result, the probe light can be diffracted into first-order and higher order directions and the diffracting power of gratings can be also enhanced at the same time by modulating the intensity Ωs of coupling fields. It is because the phase modulation is enlarged with the increase of Ωs, which makes more probe energy deflect into the high-order directions. Therefore, these gratings can be applied to the all-optical multi-channel beam splitting and switching for optical networking and communication.

Fig. 5. (color online) Diffraction intensities of zeroth-order I(θ0) (black solid line), first-order I(θ1) (blue dashed line), and second-order I(θ2) (red dashed-dotted line) versus coupling field intensity Ωs for (a) δ1 = 0.5γ, δ2 = −1.5γ; (b) δ1 = 1.5γ, δ2 = −1.5γ; (c) δ1 = 2.5γ, δ2 = −0.5γ. Other parameters are the same as those in Fig. 3.
4. Conclusions and perspectives

In this paper, by using two standing-wave fields, we propose a scheme for coherently induced grating in a mixture of two three-level atomic species with EIT and ARG systems. Three types of phase gratings are investigated in a real atomic system with 85Rb and 87Rb. By using the interaction between the two two-photon resonances, we obtain an enhanced refractive index in the medium. As a result of intensity dependent standing-wave fields, amplitude and phase are modulated periodically in space, and then a sinusoidal grating will be formed. It is shown that different diffraction patterns can be obtained by varying the two-photon detunings δ1 and δ2. The two-photon absorption in EIT system and Raman gain in ARG system contribute to the amplitude modulation, which affects the diffraction efficiency of the induced gratings. Meanwhile, we can achieve different phase-shift values by changing the intensities of two standing-wave fields. A larger phase modulation depth can help to deflect the probe energy into the higher order diffraction directions. The intensities of two standing-wave fields have an important influence on the diffraction efficiency of gratings, thereby the diffraction intensities I(θ0), I(θ1), and I(θ2) of different orders can take their maximum values at appropriate Ωs. Therefore, it is necessary to choose appropriate parameters to achieve the required diffraction and splitting of probe beam. Thus, the proposed gratings may be employed as all-optical splitter and switching in optical networking and communication. Moreover, in addition to cold atoms, EIG can be formed in thermal atoms.[19,46] Hence, our scheme can also be performed with hot rubidium atoms.[42]

Reference
[1] Harris S E 1997 Phys. Today 50 36
[2] Fleischhauer M Imamoglu A Marangos J P 2005 Rev. Mod. Phys. 77 633
[3] Dalton B J McDuff R Knight P L 1985 Opt. Acta 32 61
[4] Gray H R Whitley R M Stroud C R 1978 Opt. Lett. 3 218
[5] ZhuY F 1992 Phys. Rev. 45 R6149
[6] Fry E S Li X Nikonov D 1993 Phys. Rev. Lett. 70 3235
[7] Scully M O 1991 Phys. Rev. Lett. 67 1855
[8] Proite N A Unks B E Green J T Yavuz D D 2008 Phys. Rev. Lett. 101 147401
[9] Hau L V Harris S E Dutton Z Behroozi C H 1999 Nature 397 594
[10] Kash M M Sautenkov V A Zibrov A S Hollberg L Welch G R Lukin M D Rostovtsev Y Fry E S Scully M O 1999 Phys. Rev. Lett. 82 5229
[11] Liu C Dutton Z Behroozi C H Hau L V 2001 Nature 409 490
[12] Longdell J J Fraval E Sellars M J Manson N B 2005 Phys. Rev. Lett. 95 063601
[13] Wang H Goorskey D Xiao M 2001 Phys. Rev. Lett. 87 073601
[14] Kang H Zhu Y F 2003 Phys. Rev. Lett. 91 093601
[15] Wu Y Deng L 2004 Phys. Rev. Lett. 93 143904
[16] Huang G X Deng L Payne M G 2005 Phys. Rev. 72 016617
[17] Ling H Y Li Y Q Xiao M 1998 Phys. Rev. 57 1338
[18] Mitsunaga M Imoto N 1999 Phys. Rev. 59 4773
[19] Sheng J T Wang J Miri M A Christodoulides D N Xiao M 2015 Opt. Express 23 19777
[20] Zhai P W Su X M Gao J Y 2001 Phys. Lett. 289 27
[21] Brown A W Xiao M 2005 Opt. Lett. 30 699
[22] Zhao L Duan W Yelin S F 2010 Phys. Rev. 82 013809
[23] Wen J M Du S W Chen H Y Xiao M 2011 Appl. Phys. Lett. 98 081108
[24] Chen H X Zhang X Zhu D Y Yang C Jiang T Zheng H B Zhang Y P 2014 Phys. Rev. 90 043846
[25] Zha X W Yuan C Z Zhang Y P 2013 Laser Phys. Lett. 10 045201
[26] Zhang Y Q Wu Z K Belić M R Zheng H B Wang Z G Xiao M Zhang Y P 2015 Laser Photon. Rev. 9 331
[27] de Araujo L E E 2010 Opt. Lett. 35 977
[28] Wan R G Kou J Jiang L Jiang Y Gao J Y 2011 Phys. Rev. 83 033824
[29] Zhou F X Qi Y H Sun H Chen D J Yang J Niu Y P Gong S Q 2013 Opt. Express 21 12249
[30] Cheng G L Zhong W X Chen A X 2015 Opt. Express 23 9870
[31] Liu J B Liu N Shan C J Liu T K Li H Zheng A S Xie X T 2016 Phys. Lett. 380 2458
[32] Xiao Z H Shin S G Kim K 2010 J. Phys. B: At. Mol. Opt. Phys. 43 161004
[33] Dong Y B Guo Y H 2014 Chin. Phys. 23 074204
[34] Liu Y Wang P Peng S Y 2013 Chin. Phys. 22 104203
[35] Yu S Liao P Yang Z Y Gu W Y 2013 Acta Phys. Sin. 62 224205 in Chinese
[36] Kuang S Q Jin C S Li C 2011 Phys. Rev. 84 033831
[37] Wang L Zhou F X Hu P D Niu Y P Gong S Q 2014 J. Phys. B: At. Mol. Opt. Phys. 47 225501
[38] Wu J C Ai B Q 2015 J. Phys. B: At. Mol. Opt. Phys. 48 115504
[39] Wang L Zhou F X Guo H J Niu Y P Gong S Q 2016 Chin. Phys. 25 114205
[40] Chen Y Y Liu Z Z Wan R G 2016 Europhys. Lett. 116 64006
[41] Chen Y Y Liu Z Z Wan R G 2017 Phys. Lett. 381 1362
[42] Proite N A Unks B E Green J T Yavuz D D 2008 Phys. Rev. Lett. 101 147401
[43] O'Brien C Anisimov P M Rostovtsev Y Kocharovskaya O 2011 Phys. Rev. 84 063835
[44] Hang C Huang G X Konotop V K 2013 Phys. Rev. Lett. 110 083604
[45] Pei L Y Lu X G Bai J H Miao X X Wang R Q Wu L A Ren S W Jiao Z Y Zhu H F Fu P M Zuo Z C 2013 Phys. Rev. A. 87 063822
[46] Dong Y B Li J Y Zhou Z Y 2017 Chin. Phys. 26 014202